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We exploit millimeter wave technology to measure the reflection and transmission response of 
random dielectric media. Our samples are easily constructed from random stacks of identical, sub- 
wavelength quartz and Teflon wafers. The measurement allows us to observe the characteristic 
transmission resonances associated with localization. We show that these resonances give rise to 
enhanced attenuation even though the attenuation of homogeneous quartz and Teflon is quite low. 
We provide experimental evidence of disorder-induced slow light and superluminal group velocities, 
which, in contrast to photonic crystals, are not associated with any periodicity in the system. 
Furthermore, we observe localization even though the sample is only about four times the localization 
length, interpreting our data in terms of an effective cavity model. An algorithm for the retrieval of 
the internal parameters of random samples (localization length and average absorption rate) from 
the external measurements of the reflection and transmission coefficients is presented and applied to 
a particular random sample. The retrieved value of the absorption is in agreement with the directly 
measured value within the accuracy of the experiment. 

PACS numbers: 



I. INTRODUCTION 

Multiple scattering of electrons and photons in dis- 
ordered systems gives rise to many important physi- 
cal effects, including Anderson localization, quantum 
corrections to conductivity, localization-induced metal- 
insulator transitions, and coherent backscatteringii 2 -! 3 -! 4 -. 
Effects of disorder are most strongly exhibited in one- 
dimensional (ID) systems. The combination of localized 
resonances and bandgaps in randomly layered samples 
opens up fresh opportunities for applications and allows 
for detailed studies of slow light and other anomalous 
group velocity effects as well as the enhanced attenuation 
associated with the long path length of localized photons. 
Due to the ease of fabrication, there are many potential 
applications for ordered and disordered millimeter wave 
dielectrics, including fundamental physics studies which 
can exploit the long dwell times associated with localized 
photons. 

In this paper, we describe observations of Anderson 
localization, enhanced absorption, and slow light in mil- 
limeter waves propagating in a random dielectric. Specif- 
ically, we utilize 100-layer dielectric stacks of quartz and 
Teflon wafers randomly shuffled like a deck of cards. The 
behavior of wave propagation in this system is a mix- 
ture of band gaps in which the random dielectric be- 
comes an almost perfect reflector, resonances which are 
associated with anomalously high transmission, localized 
electric field inside the layered system, and enhanced ab- 
sorption. The block diagram of our experiment is shown 
in Fig. [TJ A key aspect of our work is that we achieve 
these disorder-induced defects in a finite system which 
is only about four times the localization length. Such 



finite size considerations are vital in quantum studies of 
disorder-induced phenomena in quantum degenerate ul- 
tracold atoms£&££. 

There are several approaches to experimental realiza- 
tion of ID systems. Recently, dielectric multilayer— ^ 
and single mode structure a 11 ' 12 ' 13 have been used to 
study localization in optical and microwave systems. 
However, millimeter waves offer several advantages. 
First, we can measure both amplitude and phase of 
the reflected and transmitted fields simultaneously. Sec- 
ond, our setup is free-space and quasi-optical, so the 
attenuation in the waveguide is minimized. Unlike 
previous work on photonic crystals and optical res- 
onators where superluminal and slow light have also been 
observe d 14 ' 15 : 16 ' 17 : 18 ' 19 , we induce these effects via disor- 
der. Ideally, a system would have sufficiently low attenu- 
ation that both the reflection and transmission response 
of a disordered system could be measured; would allow 
for amplitude and phase of the reflected and transmitted 
signal to be measured simultaneously; have large enough 
bandwidth that a number of band-gaps and resonances 
could be studied for each sample; and would allow for 
easy fabrication of samples, so that many could be inves- 
tigated quickly. All these requirements are met by our 
quasi-optical millimeter wave system. 

Our article is outlined as follows. In Sec. [Ill the exper- 
imental set-up is described. In Sec. IIIII the transmission 
and reflection data are presented. In Sec. IIV1 the effec- 
tive cavity model is derived. In Sec. [Vj the experimental 
data of Sec. IIIII are compared both to the effective cav- 
ity model and to simulations obtained from a propagator 
matrix. Finally, in Sec. IVII we conclude. 
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FIG. 1: Experimental setup. The Vector Network Ana- 
lyzer (VNA) generates centimeter waves which pass through 
a harmonic multiplier (HM), converting them to millimeter 
waves. Other components include isolators (I), a circulator 
(C), polyethylene lenses (L), harmonic detectors (HD), scalar 
(corrugated) horns (SH) and two phase-locked local oscillators 
(LO). The thick lines refer to coaxial cables, which are only 
used for the cm waves. The internal frequency synthesizer 
is locked to a Rubidium clock, allowing for long-term stabil- 
ity. Finally, the vector receiver in the VNA down-converts 
the time-varying E-fields (Tx for transmitted and Rx for re- 
flected) to a low frequency for digitization and display. The 
sample (PC) is at the focal plane of the transmitting horn. 
The thick black lines connecting the VNA and the millimeter 
wave components represent high-frequency micrwoave coaxial 
cables. Between the harmonic mixer/detectors and the scalar 
horn antennae, the millimeter waves propagate via waveguide. 



II. DESCRIPTION OF THE EXPERIMENT 



arc held tightly in the telescope tube, inevitably there 
are small air gaps between the wafers. This ultimately 
limits our ability to fit the data. 

The measurements were performed with a vector net- 
work analyzer developed by AB Millimetre. The mil- 
limeter waves are generated by a sweepable centimeter 
wave source, i.e., microwaves, in this case from 8-18 GHz. 
For work in the W-band these centimeter waves are har- 
monically multiplied by Schottky diodes, coupled into a 
waveguide, and eventually radiated into free space by a 
scalar horn. These arc also known as corrugated horns 
and have the property that when the TE01 rectangular 
waveguide mode is coupled into the circular horn, a hy- 
brid mode is established which has no cross-polarization. 
This results in the E-field being vertically polarized in a 
single-mode Gaussian beam^S. 

A polyethylene lens collimatcs the beam. The random 
dielectric stack is placed in the focal plane of the quasi- 
optical system. The transmitted field is then collected by 
an identical lens/horn combination, detected by a Schot- 
tky harmonic detector and fed to a vector receiver which 
mixes the centimeter waves down to more easily manage- 
able frequencies where the signal is digitized. Reflected 
waves are also collected by the transmitting horn and 
routed via a circulator and isolator to the vector receiver. 
The source and receiver local oscillators are phase-locked. 
A more complete description of the system is given in2i 
and22. 



Our layered dielectric samples consist of stacks of 300- 
400 fim thick wafers of fused quartz and Teflon held to- 
gether tightly in a telescope tube with retaining rings. 
These disks are stock items for Teflon and semiconduc- 
tor manufacturers and are available at relatively low cost. 
With our millimeter wave vector network analyzer, a scan 
over the entire W-band (75-110 GHz) at 5000 frequency 
points takes only a few minutes. Our setup allows for 
high signal-to-noise measurement of both reflected and 
transmitted E-fields simultaneously The electric field is 
a linearly polarized Gaussian beam whose diameter is 
much smaller than that of the wafers; the beam strikes 
the dielectric stack perpendicularly, allowing a ID treat- 
ment. 

Quartz and Teflon are low loss materials at millimeter 
wave frequencies and provide a reasonable dielectric con- 
trast. The real part of the dielectric constants are mea- 
sured to be e' = 3.80 ± .02 for quartz and e' = 2.05 ± .02 
for Teflon; these values, obtained by fitting the transmis- 
sion data for homogeneous samples with a Fabry-Perot 
model, are quite close to reference values (3.8 and 2.1, re- 
spectively^). The same fitting procedure gives loss tan- 
gents of 0.0007±0.0002 for quartz and 0.0003±0.0002 for 
Teflon. Such small loss tangents are difficult to measure 
with our setup and are less accurately specified. 

The micrometer-measured thicknesses of our individ- 
ual wafers are 0.460 ± 0.01 mm for quartz and 0.375 ± 
0.005 mm for Teflon. However, even though the wafers 



III. EXPERIMENTAL DATA 

In Fig. [21 we illustrate the transmission and reflection 
spectra for the W-band as investigated experimentally 
and theoretically for a representative sampled Plotted 
are the transmission coefficient (T) and one minus the 
reflection coefficient (R). 

The actual data collection requires careful calibration 
of the system, in the instrument response and in all phase 
shifts in the cables and waveguide. First, one sweeps over 
the band of frequencies with no sample present; this gives 
a calibration E-field Eq measured for transmission. Sec- 
ond, one sweeps with a "perfect" reflector, i.e., a polished 
metal plate, at the focal plane; this gives a calibration 
E-field Eq. Finally, with a sample in place, the sweep 
returns the amplitude and phase due to the presence of 
the sample relative to a flat instrument response. So, in 
Fig.d R = |£7££| 2 and T = |£*/£§| 2 arc the reflection 
and transmission coefficients, respectively. 

The theoretical curve in Fig. [5] is obtained by solv- 
ing Maxwell's equations numerically using the standard 
propagator matrix approach 2 ^. This allows us to com- 
pute the amplitude and phase of the reflected and trans- 
mitted fields as well as the fields everywhere inside the 
sample. As can be seen from Fig. [2] most frequencies lie 
in disorder-induced band-gaps, where the reflection coef- 
ficient is near unity. Between band-gaps there are narrow 
resonances accompanied by greatly enhanced transmis- 
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FIG. 2: Measured (top) and calculated (bottom) transmission and reflection response for a randomly layered sequence of 100 
quartz and Teflon disks. Plotted are the transmission coefficients (T, dotted curves) and one minus the reflection coefficient 
(1 — R, solid curves). 



sion. Off-resonance transmission losses, e.g., around 80 
and 94 GHz in Figure H are about 1%. 

We observe in Fig. [2] differences in both amplitude and 
frequency between the measured and calculated spectra. 
In order to understand this effect, we studied periodic 
systems (data not shown), and found that the fit was 
better. We also studied many different random samples 
and did not find any systematic amplitude and frequency 
shifts. We conclude that these shifts are likely a random 
effect of the tiny, unmeasurable air gaps between wafers 
as well as of small random variations in the thicknesses 
of the wafers. 

From the data shown in Fig. [21 the localization length 
l\ oc can be estimated. Indeed, it is well known that the 
Lyapunov exponent, 7 = l^z = — (lnT)/L is a self- 
averaging quantityi^^. That is, its fluctuations are 
small and its variance tends to zero, so that 



7^ (7) > -(1/Jloc + 1/Jabe) : 



(1) 



where L is the sample length and Z a b s is the spatial ab- 
sorption rate^l. Thus, at large L the Lyapunov exponent 
7 is close to — (1/Zioc + 1/iabs) with high probability, i.e., 
in almost all samples, or at almost all frequencies in one 
sample. We call such off-resonance samples and frequen- 
cies typical. 

In Fig. [21 we find that jL at typical frequencies is ap- 
proximately 0.01. The measured loss tangents of quartz 



and Teflon correspond to absorption lengths of approxi- 
mately 35 cm for quartz and 100 cm for Teflon. The re- 
ciprocal absorption length is 27mtan<5/A, where tan S = 
e"/e', n is the index of the homogeneous material and A 
is the free-space wavelength of the probing beam. The 
absorption lengths of both quartz and Teflon are much 
larger than the sample itself. Thus, we find for our 42 



mm sample the localization length l\ oc 



7" 



10 mm. 



Since Z a b 8 >> L, its influence on typical transmission 
frequencies is negligible, 



Ttyp ~ exp(-L/Zi oc ). 



(2) 



We reserve discussion of the atypical, or resonant case to 
Sees. HVl and IV1 



IV. THEORETICAL DESCRIPTION 

As can be seen from Fig. [51 the sum of the reflection 
and transmission coefficients drops noticeably at the res- 
onances. This means that absorption plays a significant 
role; in the presence of absorption, standard analytical 
treatments'^ fall short. Therefore, we take advantage 
of an effective cavity model recently introduced by Bliokh 
et ah 13 i 25 as follows. In the random medium of our sam- 
ple, the transmission coefficient at typical frequencies de- 
termines the localization length via Eq. © and so can 
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be measured. However, along with the typical frequen- 
cies there are resonant frequencies. These are relatively 
rare, or low probability frequencies at which resonances 
occur, and the effect of absorption is dramatic. In these 
resonances, photons become partially trapped by multi- 
ple scattering. In effect, the random multiple scattering 
gives rise to an effective cavity within the sample. The 
long dwell time of these photons is responsible for the 
enhanced attenuation. In other words, the attenuation 
increases significantly at a resonance because the energy 
absorbed at each point is proportional to the intensity 
of the field, which is exponentially large throughout the 
resonant cavity. Therefore, both the reflection and trans- 
mission coefficients decrease sharply even when the ab- 
sorption length l a bs is much larger than the total length 
L of the sample. 

One can treat any particular localized state as an ef- 
fective cavity mode. Using the measurable reflection and 
transmission coefficients on and off resonance together 
with the measurable resonance width, one can make a 
fit to determine the location and width of the effective 
cavity, as well as an effective absorption length. 

Consider a random sample as a ID resonator built of a 
cavity of length l c&v bounded by two barriers of lengths 
l\ and I2, whose transmission coefficients are 



-1,2 



exp(-Zi, 2 Aoc) , 



(3) 



with Zi oc the localization length. The parameters of the 
effective cavity will be different for each resonance inside 
the same sample. We take as our starting point the reso- 
nant transmission and reflection coefficients and resonant 
width, which can be shown to be 

4 



fires = 1 _ 



T 



8u = 



2eL 



(4) 

(5) 
(0) 



respectively, where c is the speed of light in vacuum, e 
is the average real part of the dielectric constant and 
5u> = 2irSf is the resonance width. Here 




G = 2Vi^cav/^abs, (9) 

where Z a b s is the effective absorption length such that the 
intensity I decays as exp(— x/Z a bs) in an ordinary random 
sample, and the total length of the sample is 

L = h + h + l cm . (10) 

By 



First, we determine the effective cavity length l c 
definition, 



-(£-Jcav)Ao 



(11) 



Substituting Eq. (JTTJ) and Eq. (5]) into Eq. one finds 
an implicit equation for the effective cavity length: 



lc. 



1 



,-(L-Jca T )/(2Jloc) 



(12) 



This transcendental equation is analyzed in more detail 
in Sec. El 

Second, we determine the effective absorption length 
Us- FromEqs. 0-©, 



2y/H 



v — u 



(13) 



^abs 

Then Eq. (JTT1) implies 

Us = 2V?/ cav (V U-y 1 e (£-'cav)/(2* Ioo ). (14) 

Let us solve for u, v as functions of i? rcs , T ICS alone. Equa- 
tions (gl)-© yield 



u 2 + 2uv + v 2 = 



u 2 + 2uv + v 2 = 



Auv 



1-Rr 



It follows that 



(15) 
(16) 

(17) 
(18) 



Equations (|17[) and (fT8"|) are valid only for R ICS ^ 1 and 
T res ^ 0. Substitution of Eqs. |17]), (JISJ), and JT2J) into 
Eq. (fTT]) yields the absorption length: 



_ 1 _ c l±y/IQ 

^abs — — - Z- 



r Sujy/e 1 — Rvcs - T n 



(19) 



Third, we determine the length of the first barrier of 
the cavity, h. Substituting Eqs. © and (TO]) into Eq. (J7J), 
one obtains 



(J) Eliminating u with Eq. (jT5J) and solving for li, 



1 - 1 it 1 \ 1 1 1 ^ 

'1 — ~ 'cavj — tloc m — 



(20) 



(21) 



Finally, substituting Eq. (fl"2"l) into Eq. (|2Tj) , one finds the 
length of the first barrier: 



"Zlnnln 



e Suj L 



(i± Vr* 



(22) 



Equations (jl9)) . (|T2)l . and (|22|) present the absorption 
length and the effective cavity lengths in terms of mea- 
surable quantities, namely, the resonance transmission 
and reflection coefficients T res and i? res , the resonance 
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TABLE I: Measured and retrieved parameters associated 
with 5 resonances shown in Figure [2] The subscripts on the 
frequencies refer to the peaks indicated in the plot. The lo- 
calization length is 1 cm, as obtained from the non-resonant 
transmission coefficient. 




width Slo = 2ir8f, the average dielectric constant £, and 
the localization length determined via Eq. @ from off- 
resonant transmission coefficients. In this treatment we 
have not taken advantage of the extra phase information 
that is obtained by our system. The phase provides an 
extra datum at each frequency which can be used to ob- 
tain additional cavity parameters. 

Finally, it is also convenient to invert these equations 
for the ensuing discussion to obtain 

T rcs = 4T 1 T 2 /(2Z cav /Z abs Vi+T 1 +T 2 ) 2 , (23) 
i? res = 1 - (t 2 + 2Z ca v/UsVi) T res /T 2 , (24) 
T h2 = exp [- (Z cav /2 ± Ad) /Zi oc ] , (25) 

where Ad is the shift of the cavity's center from that of 
the sample, replacing l\. 

Using data on the transmission and reflection coef- 
ficients and the peak widths for the sample shown in 
Figure [2] we have calculated internal parameters of our 
disordered system. These results are shown in Table I. 
Shown in columns 2,3 and 4 are 1 — R, T and 8f, respec- 
tively, measured for the resonances indicated in column 
1. We use an extra vertical bar between columns 4 and 5 
to emphasize the fact that the measured resonance val- 
ues lie to the left of this double-bar and the retrieved 
system parameters lie to the right. The value of T/^e 
(column 5) has been retrieved via equation [TTil The loss 
tangent tan(a) = Tc/ (2irf^/i) is presented in column 6. 
The value of the loss tangent averaged over the five res- 
onances equals 8.35 x 1CU 4 . The weighted loss tangent 
for our disordered quartz/Teflon system is 5.2 x 10~ 4 , so 
that the measured and retrieved values of the absorption 
agree to within the accuracy of the experiment. Thus 
from external measurements we have retrieved two in- 
ternal parameters of the system, the localization length, 
which relates to disorder, and the absorption length. 

Finally, we note that numerical simulations were also 
performed with the standard propagator matrix method, 
as mentioned briefly in Sec. IIII1 in order to obtain infor- 
mation about the electromagnetic field inside the sample. 
For these simulations we used the precise experimental 
ordering of quartz and Teflon wafers. 



FIG. 3: Shown is the solution to the transcendental Eq. ([26 
where A = Z C av/2Zi oc and a is given in Eq 



V. DISCUSSION 

Equation (| 1 2[) requires further analysis, as it is tran- 
scendental. We can rescale it as 



A = a cxp(A) . 



(26) 



where 



A 



r/2koc : 

c 



2 Zioc e Suj V7k 



(27) 

: exp(- J L/2Z loc ) . (28) 



Equation (|2"6"|) has a real solution in terms of a Product- 
Log for a G (-co, 1/e]: 



A = — ProductLoe 



(29) 



a special function which can be obtained numerically. 

In the physical region of a £ [0, 1/e], A — ► 1 as a — > 1/e. 
Since the localization length must be smaller than the 
cavity size, it follows that A is confined to the relatively 
narrow range of \e~^ ~ 0.3 and e _1 « 0.37. Then the 
resonance width Sto must fall off exponentially as the sys- 
tem size is increased. 

Figure 2] shows the calculated spatial distribution of 
the E-field inside the medium for the same sample as 
was used in Fig. [2J for both resonant and non-resonant 
frequencies. Figure O shows the spatial distribution of 
energy at three resonances. 

At resonance, the energy density within the localized 
mode, i.e., inside the effective cavity, can be orders of 
magnitude larger than that of the incident wave. This 
huge field enhancement has many potential applications. 
For example, at optical frequencies, this increased energy 
density has been exploited to produce a random laser- 
based on a layered (ID) medium. Off resonance, the 
dielectric multi-layer presents an almost perfect reflector; 
it is only the exponential tail of the field that penetrates 
the sample. 

Slow light is associated with strong dispersion and can 
be seen in situations involving resonant transmission, for 
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FIG. 4: Localized E field within the sample. Shown is a plot of 
the calculated E-field intensity (in dB) within the scattering 
medium. By comparing the frequencies with those in Fig. [2] 
one can see the strong localization of the field at transmission 
resonances. 
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FIG. 5: Localized modes at regions of enhanced absorption. 
E-field intensity as a function of space inside the sample at 
four resonant frequencies used in the analysis in Table I. Lo- 
calization results in a varying degrees of field enhancement 
inside the random effective cavity. The stronger the field en- 
hancement, the greater the path- length-induced absorption. 



example, with defect states in a band gap2£ or with ultra- 
cold atomic gases^i In the case of localization, slow light 
and enhanced absorption go hand in hand because they 
have the same physical origin - both of them are caused 
by multiple scattering, which increases dramatically the 
photon path and the resonance dwell time. Plotted in 
Fig.[5]is the measured frequency dependence of the phase, 
4>(f), of the signal transmitted through a random sam- 



ple. At each resonance there is a pronounced decrease of 
the inverse phase derivative, (d(j)(f)/df)~ 1 . If there were 
a well defined group, i.e., a pulsed measurement, and if 
there were a well defined local wave number, then the in- 
verse phase derivative would be proportional to the group 
velocity. For example, one could divide the transmitted 
phase by the thickness of the sample to get an average 
wave number. Such an analysis, however rough, gives 
group velocities on the order of 10% c. But our mea- 
surement doesn't have sufficient bandwidth to synthesize 
pulses small compared to the size of the sample. 

The corresponding time delay associated with a reso- 
nance can be estimated to be 



a . 1 ^cav ^ ^cav 

At oc — = — — > . 

A/ c r typ c 



(30) 



Taking Z cav of the order of twice the localization length, 
2 x 1cm, and T typ ~ 0.01, we obtain At oc Z cav /c x 
1/Tt y p = 2 cm x 100/(3 x 10 10 cm/s) « 7 ns. If we were 
speaking of well-defined pulses, such a delay would corre- 
spond to a group velocity an order of magnitude less than 
the speed of light. Although our measurements are fun- 
damentally frequency-domain, the structure of the spec- 
tral peaks of the transmission coefficient can be related 
to the conventional notion of slow light in this sense. 



VI. CONCLUSIONS 

We have developed an experimental system for study- 
ing disorder-induced wave phenomena in dielectrics at 
millimeter wave frequencies. In our system we see local- 
ization even in samples that are only four times the lo- 
calization length. The deep resonances connected with 
localization allow one to observe enhanced absorption 
and disorder-induced delays that can be associated with 
group velocities much less than the speed of light in vac- 
uum. Furthermore, off-resonance, the random dielectric 
stack becomes a near perfect reflector. This gives rise 
to a transmitted E-field phase that is nearly indepen- 
dent of frequency at the band-edge. We have interpreted 
the experimental data in terms of an effective cavity 
model, which enables us to retrieve the localization and 
absorption lengths from frequency dependent reflection 
and transmission coefficients. 

This material is based upon work supported by the Na- 
tional Science Foundation under Grants EAR-0337379, 
EAR-041292 and PHY-0547845. 
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